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Abstract. We solve the Yang- Mills equations in the framework of the McLerran- 
Venugopalan model for small times r after a collision of two nuclei. An analytic 
expansion around r = leads to explicit results for the field strength and the energy 
momentum tensor of the gluon field at early times. We then discuss constraints for the 
energy density, pressure and flow of the plasma phase that emerges after thermalization 
of the gluon field. 
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The Color Glass Condensate (CGC) model [H El El HI E] has provided valuable 
tools to understand the interaction of hadrons and nuclei at very high energies. The 
idea that the rapid growth of the gluon distribution at high energy is tamed by gluon 
recombination leads to a saturated gluon density described by a saturation scale Qg. 
Due to the high occupation numbers the gluons can be approximated by a classical field 
[1]. In high energy nuclear collisions the gluon density in each nucleus is enhanced by 
a factor ~ A^^^ and it has been argued that the CGC can describe the initial particle 
production at the Relativistic Heavy Ion Collider (RHIC) [6]. 

Here, we discuss what can be learned about the energy and momentum deposited 
in the space between the two receding nuclei immediately after the collision. This is an 
important question because there is convincing evidence that relativistic hydrodynamics 
governs the evolution of the fireball starting at rather early times tq ^ 0.5 . . . 1.0 fm/c. 
The plasma phase, given by an energy momentum tensor T^^^ = {e + p)u^u'^ — pg'^" 
emerges through a rather rapid and not yet completely understood thermalization 
process from the initial gluon field after overlap, described by an energy momentum 
tensor T^'^ . Our goal here is to constrain the energy density e, pressure p and flow 
velocity v (with = 7(1, v)) at the beginning of the plasma phase, using the CGC 
model in its version first conceived by McLerran and Venugopalan (MV) [H [2] . 

We have to solve the Yang-Mills equations [D^, F'^'^] = for a current J" created 
by infinitely thin color charge distributions pi(x^) and P2(x_l) propagating on the plus 



Early Time Evolution of High Energy Heavy Ion Collisions 



2 



and minus light cone, respectively, and overlapping at time t = 0. In light cone gauge 
the field of each nucleus before the collision, A\{-iL±) and A2{'x±) respectively, (i = 1, 2), 
is transverse. In the forward light cone, i.e. after the collision, the field has light cone 
components = ±x^A and transverse components A^. A and A^ are functions of the 
proper time r = \/t'^ — and the transverse position x^, but not of the space-time 
rapidity r/. 

For small times r, immediately after the collision, the Yang-Mills equations in the 
forward light cone can be rewritten using an expansion in powers of r, A^ = r^'A'^^y 
The resulting infinite tower of differential equations can be solved recursively to arbitrary 
order in r [7j. Explicit solutions for A^ and the field strength F^'^ up to order have 
been discussed. The lowest order in r ((9(r°)) coincides with the familiar boundary 
conditions for the field on the light cone. They imply that at the earliest time, for 
r — >• 0, strong longitudinal electric and magnetic fields between the nuclei, [HI [7] 

= zg[A\, Ai,], = -ige'^[A\, Ai,], (1) 

dominate. This was also observed, e.g., in [9l[in]. The next order in r corresponds to a 
linear build-up of transverse fields. 

One can now compute the energy momentum tensor of the gluon field. The leading 
terms 0{t') at small time are the diagonal elements 

e = = T/^ = -Tf (2) 

{i = 1, 2). Here e = {E^ + Bl)/2 denotes the energy density at r = 0. It turns out that 
e in the McLerran-Venugopalan model suffers from a ultraviolet divergence [3 [9]. In 
lattice calculations this divergence is regularized by the finite lattice spacing [9]. We can 
argue that this divergence comes from the fact that the ultraviolet sector of the initial 
particle production with large momenta Qs is not perfectly described in the MV 

model which works best for ~ Qs- Instead, we can choose a cutoff Qq and replace 
the classical field above Qq with a quantum but weak coupling description known to be 
well-behaved in this regime, i.e. perturbative QCD. 

At order 0{t^) the components describing transverse fiow of energy receive their 
leading contributions 

Tf°* = cosh r], Tf = sinh r/, (3) 

with the transverse fiow vector 

y' = ~y'{El + B';) (4) 

{i = 1,2). We will see below that this translates directly to the existence of transverse 
fiow in the plasma phase immediately after thermalization. 

Further corrections to the components of T/"^ up to order 0{t^) can be easily 
computed. This gives a rather accurate picture of the energy momentum tensor at times 
r << 1/Qs- At T ^ 1/Qs the expansion might fail, however the asymptotic behavior 
for r — i> oo is known from a weak coupling expansion [2] and an interpolation can be 
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used for intermediate values of r. This reproduces the time dependence of numerical 
solutions of the Yang-Mills equations, see e.g. [9]. 

In any case, the classical field approximation is not expected to hold for very large 
times anyway. Instead, the field is expected to decay on a time scale ~ maybe 
through instabilities in the color fields pTl] or particle production in the background 
gluon field [12]. This is widely assumed, though not proved, to lead to a thermalized 
plasma within a rather short time. 

One can, however, use energy and momentum conservation to estimate the energy 
momentum tensor of the plasma. Suppose the thermalization is sufficiently rapid around 
a time tq so that the system can be approximately described by T^^'^ for t < tq, and by 
T^j'^ for T > Tq. Then one can show that 



= , V = tanhr) (6) 

cosh rj e + p 

{i = 1, 2) using the first two terms in the r expansion. These four equations relate 
the five parameters e, p and v in the plasma phase with components of the energy 
momentum tensor of the field. The system of equations can be closed by providing an 
additional constraint, e.g. an equation of state. This result can be used as an initial 
condition for a hydrodynamic evolution of the system. 

Let us briefly discuss this result. First we notice that we recover boost-invariance 
for the longitudinal flow velocity v^. Secondly, as mentioned before, there is initial 
radial flow f * in the plasma phase which is directly proportional to the energy flow z/* 
of the field. For collisions with finite impact parameter this also implies the existence 
of elliptic flow. 

The results presented so far are functions of the single nucleus fields A\ 2 and have 
to be evaluated by putting the correct expressions for those fields. They have been 
discussed in the literature, see e.g. [3j. In |7] estimates were presented in an abelianized 
approximation in which the non-abelian effects were mimicked by color screening at a 
typical distance Rc = l/Qs- To be precise, the screening radius Rc at a given point in 
the transverse plane was chosen to to be R~^ = iasCr/S where a is the number density 
of color charges in the nucleus [7j. This approximation reproduces the gluon two-point 
function {A\x)A'^{y)) [3] quite well. Thus we expect this approximation to give a good 
estimate of the energy momentum tensor, since it has been shown that e depends on 
A\ 2 solely through this two-point function [9] . Using this scheme one finds very simple 
formulae for the center of two very large nuclei colliding head-on, 

e =^-^a,a,\n{l + ce), (7) 

u^ = -r^V\a,a,)\n{l + ce), (8) 
where c ~ 0.42 is a numerical constant and ( = RcQo- 0^1,2 are the number densities of 
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charges in the two nuclei, which lead to the SU{3) color densities pi,2 [Z]- 

Fig. 1 shows the dependence of the results on the (unphysical) cutoff Qq. The 
(physical) saturation scale Qs = is independent of Qq as it should be, while the 
initial energy density e (eo in the plot) grows with Qq. For a rather reasonable value 
Qq = 2.5 GeV the estimated initial energy density is roughly e = 260 GeV/fm^. 
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